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ABSTRACT 

The vaporization of a droplet, interacting with its neighbors in a non- 
dilute spray environment is examined as well as a vaporization scaling law 
established on the basis of a recently developed theory of renormalized 
droplet. 1 The interacting droplet consists of a centrally located droplet and 
its vapor bubble which is surrounded by a cloud of droplets. The distribution 
of the droplets and the size of cloud are characterized by a pair-distribution 
function. The vaporization of a droplet is retarded by the collective thermal 
quenching, vapor concentration accumulated in outer sphere, and by the limited 
percolative passages for mass, momentum and energy fluxes. The retardation is 
scaled by the local collective interaction parameters; group combustion number 
of renormalized droplet, droplet spacing, renormalization number and the local 
ambient conditions. The numerical results of a selected case study reveal 
that the vaporization correction factor falls from unity monotonically as the 
group combustion number increases, and saturation is likely to occur when the 
group combustion number reaches 35-40 with interdroplet spacing of 7.5 
diameters and the environment temperature of 500 K. The scaling law suggests 
that dense sprays can be classified into: (D a "Diffusively Dense" cloud 

characterized by uniform thermal quenching in the cloud, (2) a "Stratified 
Dense" cloud characterized by a radial stratification in temperature by the 
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differential thermal quenching of the cloud, or (3) a "Sharply Dense" cloud 
marked by fine structure in the quasi-droplet cloud and the corresponding 
variation in the correction factor due to the variation in the topological 
structure of the cloud characterized by pair-distribution function of quasi- 
droplets . 

NOMENCLATURE 

D Mass diffusion coefficient 

g Pair - distribution function 

G Group combustion number 

L Latent Heat of vaporization 

m Vaporization rate 

n Number density 

q Heat of combustion 

r Radial coordinate 

R Radius 

s Droplet separation 

T Temperature 

u Velocity 

W Molecular weight 

a Schvab-Zeldovich Variable 

6 Renormalization number n /n 

ts co 

Y f L/q 

e F - ( Vf )_1 

6 r £ (n)/r^(o) 

v Stoichiometric coefficient 

n r/r 

l 

p Density 
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Subscript 


c Canonical bubble 

co Canonical bubble of test droplet 

F Fuel 

T Temperature 

i Liquid drop surface 

ts Transition sphere surface 
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1. Introduction 


Increasing theoretical findings ^ ^ and experimental evidence 14-18 
attest the widely held belief that the short-range collective interaction 

H-13 

among the neighboring droplets and the long-range interaction l > 7 > 9 , 
with the droplets at distance, on a hydrodynamic scale, have profound 
impact on the state of a droplet, i.e. the states of saturation, vaporization, 
ignition, combustion and extinction, as well as the droplet interfacial 
process rates in non-dilute cloud or spray environments. These collective 
interactions, produced by local hydrodynamic and transport processes in a 
complex topological environment, serve to control the percolative passage for 
dispersing mass, momentum energy fluxes and the effective interfacial area for 
the property exchange processes. These interactions result in the collective 
thermal quenching, the accumulation of vaporizing species and the tendency for 
stagnating the microscale local Stefan flow and mean flow through dynamic 
equilibration between the two phases. 

Review of the current theories of collective interactions including: 

Group Vaporization and Combustion 7 > 10 (GVC), Discrete Droplet Model (DDM) 

2 3 5 

’ ’ and Droplet In Bubble (DIB) 4 > 8 reveals two major theoretical 

deficiencies in the theory of short-range interaction and droplet rate 
processes at this juncture. These are: (1) the lack of the fundamental 

concepts and the mechanisms interlinking small scale discrete droplet 

processes with a large scale quasi-continuum flow of a non-dilute cloud or a 

spray, and (2) the incomplete understanding of complex interfacial processes, 
finite rate reaction, turbulence and transport processes that occur in the 
vicinity of each droplet. The first issue compels the current research of 
non-dilute sprays to proceed with two rival approaches; GVC, that primarily 
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deals with the long-range interaction, thus neglecting small scale resolution, 
and DDM, which focuses attention on the detailed behavior of discrete droplets 
in well-ordered droplet assemblies, thereby requiring a laborious compu- 
ational procedure involving a large number of droplets encountered in 
practical sprays. 

The alternative approach, adopting a continuous spray model supplemented 
by realistic droplet laws, appears the compromised mean of the prediction of 
non— dilute spray. A recent study of Tishkoff ^ on the numerically correlated 
vaporization correction factor derived from the DIB model demonstrates the 
viability of such combined approach to complement modern non-dilute spray 
calculations formulated on either Eulerian— Eulerian or Euler ian— Lagrangian 
framework. It must be mentioned, however, that an attempt to deduce an 
improved droplet vaporization law from the results of existing DDM met with 
difficulties due to the geniune lack of a self-consistent criterion of 
"droplet environment" whereupon the gas properties such as the temperature and 
concentration of vapor species, are inserted in the droplet law for the 
determination of the vaporization rates. A universal theory of short-range 
interaction and a set of comprehensive laws of droplet rate processes remain 
as the major unsolved issues in the contemporary theory of non-dilute sprays 
which has been the central theme of research 13 conducted. The objectives 
of this paper are to present the basic concepts, theoretical approaches and 
the results of Renormalized Droplet (RND) theory to establish vaporization 
laws for droplets in a stationary non-dilute cloud environment. 

The paper begins with the descriptions of droplet models and theories, in 
section 2, to clarify basic concepts and definitions of "droplet," adopted in 
modern spray theory. The structure, model and mathematical analysis of RND 
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are described in Section 3. The results of numerically predicted RND 

structure, vaporization correction factor and their dependence on collective 
interactive parameters are presented in Section 4. 

2, Droplet Models: Topological Properties and Modeling 

2. 1 Elemental Model 

Three fundamental topological properties that play basic roles in 
interaction phenomena between the droplet under investigation, termed the 
,! test droplet", and its neighbors, i.e. "field droplets," are; (1) 

graininess or discreteness" of a droplet, (2) "covolume" or "evacuated 
volume" of a droplet, and (3) "localization" of a droplet relative to a 
reference point. Basic droplet models which attempt to capture properties of 
the droplet to a desired level of accuracy, are classified into three types 
according to the level of sophistication in the characterization of the 
topological features, as summarized in Table 1. First, the "natural droplet 
model" which characterizes a droplet by its true geometrical shape, size and 
location, simulates all the topological properties described above. 
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Table 1 
Droplet Models 


Droplet 

Properties 

Natural 
Drop Model 

Quasi -Drop let Mode 

;i 

Non-Unif ormly 
Smeared Drop 

Uniformly 
Sheared Drop 

Point Source 
Model 

Graininess 

Sharp* 

Diffuse** 

No 

Singular 

Covolume 

Sharp 

No 

No 

No 

Lo calization 

Sharp 

Diffuse^ 

No 

Singular 


* Sharp properties caused by phase discontinuity 
**Diffuse properties featured by the absence of phase discontinuity 


The natural droplet model provides a sharp droplet configuration required 
for the predictions of microscale flow structure, and results from a 
combination of convective and Stefan flows, and interfacial process rates of 
the test droplet. The solution is given by solving Dirichelt boundary value 
problems associated with conservation equations and a well defined set of 
physical boundaries formed by the interface of droplets. 

The second type of model, i.e. "quasi-smeared droplet," describes a 
droplet by a medium spreading through the space. The droplet thus coexists 
with the host medium without apparent phase discontinuity. The models in this 
category are further divided into "uniformly" and "non-unif orraly" smeared 
droplets, Fig. 1. In the latter model,, the location and spatial extension of 
droplets are depicted by a joint probability distribution function. The 
"graininess" and "localization" of droplets are therefore partially 
characterized by the joint probability function rather than the phase 
discont inuities as in the case of the natural drop model. 

"The uniformly smeared droplet model" which has been adopted in the 
majority of two-phase flow and spray theories is a special case with a uniform 
distribution of the joint probability with its numerical value equals to unity 
throughout the space. In this model, "graininess" and "localization," are 
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absent. Because of the lack of the basic topological properties, these quasi- 
droplet models are unable to be used for the determination of the droplet 
bound flow field and droplet laws. However, when appropriate droplet laws are 
provided to describe the interfacial processes rates, the flow structure of a 
farfield can be predicted with an acceptable accuracy by a quasi-two phase 
approach that is simpler than the Dirichelt boundary value problems. 

The third model is a "point source model" which uses a mathematical 
singularity to characterize the "location" and a "singular graininess" of a 
droplet in a dilute spray with a large spacing, i.e., the ratio of spacing to 
the droplet size is much greater than unity. Like a "quasi-droplet model", 
the point source model" fails to provide a satisfactory near-field structure 
required for the prediction of the rates of interfacial processes, but the 
model offers a simple mathematical theory for the prediction of far field 
solutions. 

2.1.2 Composite Hodel 

Although a specific elemental model has been frequently used in the 

analysis of single droplet, or many-droplet problems (e.g. natural droplet 
model in DDM and DIB, uniformly smeared droplet in GVC), the potential 
advantages of the simultaneous use of two or more than two elemental models 
has not been fully exploited. 

A composite modeling technique permits a strategic selection of a desired 

elemental model in a certain selected region and an alternative model in 

another region to enhance the modeling flexibility and the simplicity in 
analysis. The choice of an elemental model is determined by the type of the 
flow field data and accuracy required in each region. For example in RND 
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theory, the test droplet is modeled by a "natural droplet model" to facilitate 
the determination of the near flow field, which depends on the detailed 
geometry of the test droplet, and the flow disturbance created by the 
neighboring field droplets# On the other hand, "field droplets" are modeled 
by "non-unif ormly smeared droplet" which are distributed in the neighborhood 
of the test droplet# The distribution of these quasi— droplets is described by 
a joint probability function to simulate localization and graininess of the 
quasi-droplets. This composite model provides a simple, self-consistent and 
useful analytical method of treating the interfacial phenomena of a test 
droplet interacting with its neighbors# The details of the application of 

composite model in RND is described in the remaining part of this section and 
the next section. 

2*2 Theories of Droplet: Canonical and Renormalized Representation of a 

Droplet Under Short Range Interaction* 

"The Theory of Droplet" concerns itself with interfacial phenomena and 
the process rates of a test droplet in isolation or under the influence of the 
collective interaction of field droplets# A representative analytical 

approach for the former is the single droplet theory, and for the latter case, 
the approaches are DIB and RND models# Since the latter problem concerns 
collective interaction between a test droplet and field droplets, the nature 
of the problems and the theoretical procedures are, in general, similar to 
that of "many droplet problems". However, an important distinction between 
the two theories is that the "droplet theory" aims to establish the rate of an 
interfacial process of a test droplet on the framework of "the minimum^sized 
many— droplet system" that includes the smallest number, N m j_ n , of droplets 
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i.e., a test droplet and N in -l field droplets confined in, "the minimum sized 
region", i.e. region. The requirements of "the minimum number of 
droplets" in "the minimum region" are imposed in RND theory so that the 
droplet laws, deduced from such a minimum droplet system, can be explicitly 
formulated in terms of the self “consistently defined "local properties of a 
spray flow field. This basic feature of "local representation," is essential 
in modern spray calculations that use local droplet laws. In contrast to the 
theory of renormalized droplet, the "cluster model^’ treats a droplet 
system with an arbitrary number of droplets in a finite region. Thus, the 
model is appropriate for dealing with a partial or complete domain of a spray 
by an approach that is different from the conventional spray theory. 

2.2.1 Canonical Droplet Theory^: CDT 

The DIB^ which consists of a natural droplet and its Wigner-Seitz bubble, 
has two fundamental properties. Firstly, a high degree of symmetry preserved 
in the DIB's periodical droplet assembly inhibits the transports of mass, 
momentum and energy transports among the neighbor droplets. Thus, DIB 
portrays a thermo — chemically closed adiabatic system with respect to neighbor 
droplets. Secondly, the strength of the interaction between the droplet and 
its bubble is determined by the size of the bubble which has the radius equals 
to half of the droplet spacing. When the droplet separation becomes 
infinitely large, the predicted vaporization rate approaches to that of an 
isolated droplet, as expected. 

Because of these two unique features of the DIB which may be regarded as 
a reference model of a short-range interaction, the theory will be termed 
"canonical droplet theory" (CDT). The theory provides a basis of RND model 
described in the following subsection. 
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2.2.2 Renormalized Droplet Theory 

Since the canonical theory treats the interfacial processes of a droplet 
in a hypothetical, well ordered droplet lattice, the validity of the canonical 
droplet law is questionable when applied to practical sprays with disordered 
droplet distributions. Needless to say, the adiabaticity and closure of the 
test droplet relative to its neighbor also break down in practical sprays. 
Another theoretical deficiency of DIB is the lack of the geometrical 
compatibility of the "environment" with that of droplet environment in spray 
theory. In CDT, the edge of a bubble has been adopted as the environment, 

though such choice is not necessarily the unique alternative, and the gas 

temperature or species concentration is used for the determination of the 
vaporization rate and droplet transient heating rate. However, in a spray 

calculation, the local average gas phase temperature of multi-phase flow is 

used for the determination of droplet process rates. Clearly, the environment 
of a canonical droplet does not coincide with that of a spray. Such 
incompatibility prohibits the encoding of the droplet law, derived by DIB, 
into a spray calculation unless an explicit link between two environmental 
properties, i.e., DIB vs smeared droplet model adopted in a conventional 
multi-phase continuum spray is provided. 

The RND theory provides a rigorous theoretical procedure that removes the 
major shortcomings of the canonical theory and provides encodable droplet laws 
by the applications of (1) composite droplet modeling technique and (2) the 
minimum sized many-droplet system, as described in the following section. 
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3. Renormalized Droplet (RND) 

3. 1 Structure of Model 

The theory of renormalization^ portrays a test droplet interacting with 

its neighboring field droplets by a composite "dressed or clothed" droplet 

model shown in Fig. 2. Two principal structural elements of RND are (1) a 

Drop let -In-Bubble (DIB) and (2) a cloud of non-unif ormly smeared droplets 

which functions as an external clothing for the DIB. The distribution of the 

quasi-droplet is described by the pair-distribution function representing the 

joint probabiltiy of finding a test droplet and its neighbor at a separation 

s. The pair-distribution function vanishes in the immediate vincinity of the 

bubble, representing the evacuated volume effect of droplets, and then 

increases rather rapidly to a value greater than unity at the radial distance 

comparable to a mean droplet spacing. This first high droplet density region 

populated by the nearest neighbor droplets is termed the "first coordination 

shell." The population peaks of the "second and higher order coordination 

shells" lose their sharpness as they merge with one another and are ultimately 

lost in a continuous environment where the pair distribution function 

approaches to unity. The size of the transition sphere, R ts , defined as the 

radial location where the pair distribution function is 0.99, depends on the 

droplet size, spacing and arrangement. The ratio of the size of the 

transition sphere, R , to the characteristic hydrodynamic scale, L, is 

t s 

typically much smaller than unity, and thus the sphere degenerates into a 

point in the limit when R /L vanishes. Accordingly, the average properties 

ts 

of the gas over the transition sphere approach to the local properties of the 
continuum flow as the limit R /L goes to zero. This "correspondence 
hypothesis" and the interlinking transition sphere constitute the two key 
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factors for the determination of the droplet process rate by the local gas 
properties of the continuum flow field. In the present analysis, RND is 

assumed to be spherically symmetric, and non-reacting. Transport properties 
are constant and the Lewis number is unity. 

The criterion of a quasi-steady state is a much more complex issue than 
that of a single droplet because of the multi-scale and multi-time phenomena 
linked with mass and energy transport in typical RND. In general, the quasi- 
steady state assumption is valid when (1) the largest characteristic diffusion 
and conduction times associated with the transition sphere is much smaller 

than the life time of any droplet in the cloud, (2) no droplet in the 

transition sphere is in the state of transient heating, and (3) no gas phase 
region is experiencing an initial or a terminal transient process. The 
effects of transient processes, and the validity or the limitation of the 

quasi-steady theory in a practical spray are discussed in Section 5. Detailed 
time-dependent analysis of transient processes in RND, which will be presented 
as Part II In the future, reveal that RND is found in some dense sprays to 
exhibit only a brief or finite period of quasi-equilibrium state. 

Additionally, RND exhibits dynamic saturation in time scales comparable to the 
characteristic diffusion time in the canonical bubble when the local group 
combustion number exceeds a critical value. 


3.2 Mathematical Analysis 

Non-dimensional equations governing RND are formulated separately for the 
inner DIB region and an external quasi-droplet cloud, respectively, as 


follows. 


,diK 


1 d , 2^1, n 

— sr (n sr’ _ 0 


1< n < n 

— — CO 


(l) 
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and 
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7“-) = G gpB. 
dn si 
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pu 
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CO 


<_ n <. n 


ts 


( 2 ) 

( 3 ) 



where n = r/r^(0), p and v are non-diraensionalized by the gas density, and 

3 

velocity on the droplet surface, respectively, Gs = 4imr , g is a pair 

9,0 

distribution function, \i is the vaporization shape factor defined by (12), 
iJj . are the fluxes of various properties; mass for i = m, fuel vapor for i=F 
and thermal energy for i=T, and A refers to the properties pertaining to an 
outer region. The definitions of properties cu and constants B^ are 

summarized in Table 2. 


Table 2 

Schvab-Zeldovichs Variables and Constants 


) 


■ 

a . 
i 

a . 

l 

B. 

l 

M 

-<V u fV 

iVw 

e F = - (W F V 1 

m 

r r c p dT/q 

b 

C TdT/q 
b 

Y F = L /q 


The system of Eqs. (l)-(4) is integrated by a repeated quadrature. The 
constants of integration are determined from the conditions of the 
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impermeability of the gaseous species on the droplet surface and the balance 


of the heat conducted to the droplet with the latent heat of the 
vaporization. The integration gives the following vaporization laws 


m(0) = 


4TrpDr^(o) 

"T^7 


JLn ( 1 + 


Tco 


-) 


(5) 


co 


or 


m(o) = 


rirpDr^ (o) 


in (1 + 


Fco - 
a F i 


x Fl 


- e. 


1 - C 


co 


where 


( 6 ) 


-1 

5 = rt 

^co co 

These laws, which agree with those obtained by Tishkoff^, are valid for 

any value of t , provided the bubble contains no droplet other than the test 
co 

droplet. This standard law will be termed "canonical vaporization law, the 
spherical bubble surface will be identified by Canonical environment’* and the 


bubble temperature 

by "canonical 

temperature • " 

In 

contrast to 

what 

is 

described above, 

an alternative 

law in which 

the 

vaporization 

rate 

is 


determined by the gas temperature on the surface of the transition sphere (see 
Fig. 2) will be called "renormalized vaporization law,*' the surface of 
transition sphere is ’’renormalized environment** and temperature on the surface 
is the "renormalized temperature" which is numerically equal to the local gas 
temperature of the continuum flow field. 
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Renormalized Droplet I^w 


Since the continuum theory of sprays predicts the local gas temperature 

but not the canonical temperature, the canonical vaporization laws (5) and (6) 

cannot be used in spray calculation as previously described. The alternative 

proposed ia this paper is to use a renormalized vaporization law. The 

renormalized law is described in this subsection. By adopting a mathematical 

procedure involving appropriate linear combinations of Eq. (2) governing 

^ ^ and iJj , one obtains two homogeneous equations governing 

m F T 

a + y and - e . These two equations are integrated and joined with the 
T F F F 

inner solutions on the surface of the canonical bubble. The resulting 
solutions are given by 


a T ^ F a F e F 

a Tts Y F a Fts e F 

, n ts - ? 2 2 

= exp { / — r ( 1 + G f gy£ “df; )dn '} (7) 

O' T] ^ S 


In order to obtain the renormalization laws, one first determines 

a and by replacing the n appearing in the lower limit of the 

Tco Fco 

integration of Eq. (7) by n • Subsequently, by substituting the resulting 

co 

expressions of and a„ in Eq. (5) and (6) one arrives at the following 

Tco Fco 

renormalized vaporization laws: 


Tt ci 

m(o) = 47rpDr (o)C ln(l + — ™) 
£ v y r 


( 8 ) 


a Fts a F£ 

m(o) = 4-npDr^ (o )C v In (1 +- ) 

a FP, “ e F 


(9) 
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correction factor calculated as follows 


where C is the vaporization 
v 



in which 


( 10 ) 


K(n|c) = n 2 C 2 g(0 


a Tc (n ) 

iCn) 4 ,pDr t ln|l *TT 1 


m(o ) m(o) [1“5 C Cn ) ] 


( 12 ) 


where a F is given by 

/S A 

a F ts + a Tts + e F (13) 

a F£ 

a Tts + Y F 

The vaporization laws (8) and (9) remain incomplete until the distribution of 
p is provided. The determination of p requires the knowledge of the 
canonical temperature a Tc (n) of a quasi-droplet located at n, see Eq. (12). A 
theoretical procedure of the determination of « Tc (n) in terms of the local gas 
phase temperature, a^,(n), is described in the following. 


Mean Canonical Bubble Temperature 
of a Field Droplet 

In the analytical estimate of a mean canonical temperature one assumes 
that each field droplet in the transition sphere of RND has a canonical 
structure, i.e., the temperature and concentration profiles are those given by 
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canonical model. Consider a canonical field droplet located at the field 
point, n , measured from the center of the test droplet. Let y be the radial 
location of an arbitrary point within the bubble of the field droplet measured 
from the center of the field droplet. Then, according to the result of the 

DIB model, the temperature distribution in the bubble of the canonical field 
droplet is given by 


« T (n,y) = y f {exp l~~ (1 - i)-l]} 

(n) y 


(14) 


where 


a = 4 it p D r £ (n) 

Additionally, the canonical temperature, cx^Cn) is obtained by replacing 
y in Eq. (14) by y c ; i.e., the radius of the bubble. The result is 

a Tc(n )” a T ( "') , c ) * V 6 *” <» - f-» - 1> os) 

c 

Since the local gas temperature a T (n) equals to the volume averaged mean gas 
temperature in the bubble of a field droplet, one writes 

- 1 y c 

a T (n) = ~ f 1 a T (n,y)y dy 

C (16) 
where V c is the non-dimensional volume of the bubble given by 
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inroduces the 


v = 4 (y c -1) 

c 3 c 


In order to express a Tc (n) in terns of a T (ri)> one 


correlation function X as 


A(n,y ) = a Tc (n,y c )/a T (n) 


(17) 


By comparing Eqs. (16) and (17), one identifies A(ri,y c > as follows 


X(ri,y c ) = 


{expl ^T (l ^" ) ' 11 ' 1} (y ^ 


3 / y i C{exp [ 7(4 t ( 1 -7 )] ' 1} ydy 


(18) 


Thus, the non-dimensional vaporization rate, u , of a field droplet is given by 


. , . 4upDr 
_ m(n ) _ * 




(19) 


m(o ) m(o ) l“C c ( r >) 


where z, ^ (n)=^ cQ 0 (n )g(n ) 


1/3 


(20) 


substituting (19) into Eq. (2), for one obtains the following 


By 

equation. 


ijt+iil-c 8g 0 

dn 2 n dn s 1 ^ (n) 


( 21 ) 


where 


83 


( 22 ) 


d *m a T (n) 

T-- = £ n (1+ — ) 

d n Y p 


In the ensuing analysis, X will be assumed to be unity. This corresponds 
to droplet vaporization with a high transfer number. The general case, 

A* l, can be solved by an iterative analysis with a guessed value of o(n) to 
calculate approximate <f> from Eq. (21). The iteration will continue until the 
vaporization rate prediced from the iterative solution converages to the 
guessed value of a. It is expected that the temperature rise in the radial 
direction is faster when A>1. Details of the general case will be reported in 
the future. 

The last step required for determining p which appears in the scaling law 

is to determine 4> as a function of q. With A=l, one can show that the 

expression for * which satisfies at n = n cQ and at n =n ts is given 

by a linear combination of two homogeneous solutions W and W , that satisfy 

dw dV, 

the canonical boundary conditions; W,(n ) = 1, - (n )=0 and ~(r\ ) = 1 

1 CO dn CO J dn v 'rr> » 


w C n ) 

* M ■ ho {u i (n) + wTT-> 

CO 2 ts) 


(23) 


where 4> cq , the characteristic value, is determined in terms of * as 
follows. 

By equating the vaporization rate in the canonical form (5) and 
renormalized form (3) and by using the definition of <Ji given by Eq. (22), one 
obtains 


4irpDr (o) 

m(o) = <j> 

1-C _ T co 


AirpDr^ (o) 


K ts 


CO 


n ts n 

1 ~ z ts + G s ^ K (h|c)h(e)dEdn 

CO CO 


(24) 
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where y (C ) -0(C) 


( 25 ) 


6(0 = r t U)l r iQ 


On substituting the expression of <t>(h) given by Eq. (23) into RHS of Eq. 
(25) and by inserting the resulting expression in the p-term appears in the 


denominator of the Eq. (24), one obtains 


r--r a - A(G s'"co-\s> 

ts ts 


where 




^co W 2 ( ° 


CO CO 


r ' ts 


i-c t .« f ts / Kt-ilOoCOr^y !«!<«> “2 


V W n (e)]dgdTi 


n n 

CO CO 


Finally, the scaling factor Eq. (10) can be expressed in terms of 


and ns follows 


{l - C t8 + G s L tS /" Wn|O0(5) [»!<«) + BH 2 < £ >]de d n} 

' CO CO 
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4 Numerical Results 

A numerical analysis will examine the structure and scaling laws of RND 
and their dependence on the principal collective interaction parameters of a 
stationary cloud of n-octane droplets with the following fuel properties; 
p£ -707 Kg/m T b = 398.7 k, L = 71.7 Kcal/Kg, and W p = 114.14 Kg/kg-mole. 

(1) Pair-Distribution Function and Canonical Bubble 

In the absence of experimental data, a pair-distribution function is 
constructed on the basis of the geometrical distribution of molecules in a 

dense liquid. The following two parameter function is adopted for the 
numerical calculation. 


g(h) 1 + a exp (-bn) cos (27rr, co ), 2n co < n < n t8 (29) 

where a and b are constants to be determined from the experimental data. In 
the present analysis a = 1.8, b = 0.65 are chosen. The resulting distribution 
patterns are shown in Fig. 3. The corresponding signatures of the inverse of 

the radius of canonical bubbles with n cQ = 5, 10, and 15 are shown in Fig. 

4. 

(2) Vaporization Shape Factor 

A pronounced increase of v in the radial direction is observed for a 
smaller value of droplet spacing, i.e. n cQ = 7.5 in Fig. 5. The ratio of the 
vaporization rate of the droplet located at the first coordination shell, for 
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the case of n = 7.5, is approximately 5% of the corresponding value of RND 
co 

with n = 15. This trend of a higher increasing rate of u at smaller 
co 

canonical bubble radii, or droplet spacing, is a common feature for small 
droplet spacing. This is confirmed for the cases characterized by pair- 
distribution functions with different values of a and b. The results suggest 
that "y -stratification" is a unique feature of RND in non-dilute clouds or 
spray. 

(3) Temperature Distribution 

High p- stratification at smaller droplet spacing, as illustrated in Fig. 

5 is due to the steep radial temperature gradients in the transition sphere, 
shown in Fig. (6). Indeed, the comparison of "u - T stratifications" suggest 
that (1) the rapid vaporization in outer layer of the cloud collectively 
quench the environment and thereby reduce the vaporization of the test droplet 
and (2) the increase in n ts at a fixed droplet spacing tends to reduce the 
inward heat transfer rate and thus suppresses the vaporization of the test 

droplet. 

(4) Vaporization Rate - Correction Factor 

The correction factors of RNDs for three selected values of n cQ , Fig. 7, 
are found to decrease monotonically as the group combustion number of RND 
increases. Saturation is projected to occur when G RN =30 ~ 40 with n co = 7 * 5 * 
While the group combustion number is a primary factor controlling the 
magnitude of C v> Fig. 7 shows a small variation in the correction factor for 
two RND's which have the same group combustion number but different 
renormalization number, 8 = /h c0 » 
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4 


Discussion 


Scaling Law with Linear Stratification Model 

In order to facilitate the practical application of the scaling law, the cor 
rection factor is integrated by adopting the following functional form of g p 


gp = guKjT) 


where is a stratification coefficient, and gp is the mean value of the 
weighted vaporization shape factor. The correction factor predicted by the 
linear stratification model described above, is given by 


C V {1 C ts +G RN [(1 ' 6 )(1+ 1 " p 2 ) + 1 K l\s (1_ ?>< 1+ | + "T " 'V ]} (30) 
where G RN = j G^iT. 


SPRAY CLASSIFICATION 

Numerical assessment of the scaling law suggests the following structural 
classification of non—dilute sprays# 

1* Diffusively Dense Cloud 

In a moderately dense cloud, RND is expected to have a large transition 
sphere that has no p- stratification. The reduction in the vaporization is 
attributed to uniform thermal quenching. The correction factor given below is 
described by the group combustion alone: 
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2. Densely Stratified Cloud 

This cloud is featured with a strong stratification in a transition 
sphere that causes an intense collective quenching and a reduction in the 
vaporization rate of the test droplet. The renormalization number ti Cg /n c0 is 
larger than unity so that the correction factor is given by 


C 

v 


1 - 5 


ts 


g rn (1 + I K 1\ S >!' 1 


(32) 


3. Sharply Dense Cloud - Fine Structure 

When the coordination shells contain the largest possible number of the 

droplet so that the renormalization number is not excessively large compared 

with unity, the correction factor depends on all the collective parameters: 

G , n /ri and n • Two sharply dense clouds with the same group combustion 
RN ts co ts 

number will exhibit structural variation when the renormalization number is 
different. 

5 Conclusion 

The present theory portrays a droplet interacting with its neighbors by a 
minimum sized "dressed" droplet structured with a centrally located droplet 
and its bubble, e.g. Droplet -in-Bubble , surrounded by a quasi -droplet cloud 
spreading through a transition sphere. The distribution of quasi —droplets and 
the size of the cloud are described by a pair distribution function. The 

introduction of the transition sphere and the correspondence hypothesis 
postulating the equivalence between the average gas properties on the 
transition sphere with the gas properties of the continuum flow field, as the 
limit R ts goes to zero, constitutes the fundamental link between the 
discrete droplet dominated region with that of the surrounding continuum 
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flow* 


This feature provides the basis of the theory of short-range 


interaction and droplet rate processes* 

The vaporization scaling law reveals that the droplet vaporization is 
retarded primarily by the collective thermal quenching and the formation of 
high vapor concent rat ion in the transition sphere. The retardation is scaled 
by combined topological-therraochemical parameters: the group combustion 

number, renormalization number and droplet spacing. The scaling laws suggest 
that non-dilute clouds or sprays can be classified into: (1) "Diffusively 

Dense" clouds in which the vaporization reduction is attributed to the effects 
of uniform thermal quencning that can be scaled by the group combustion number 
alone; (2) "Densely Stratified" clouds, with marked radial stratification due 
to the large temperature gradient which exists in the transition sphere; and 
(3) "Sharply Dense" clouds, with a relatively large canonical bubble, and 
radial stratification so that the scaling law depends on G R ^, r\^_ and 8. 
Selected case studies show that the correction factor falls monotonically as 
the group combustion number increases and the saturation is expected to occur 
when the G RN is of the order of 30 - 40, with the droplet spacing of 
approximately 7.5 times of the droplet size and environment temperature of 
500 K. 

Finally, the validity of quasi-steady (Q-S) approximation is examined by 

comparing the relative magnitude of the characteristic time, t , of 

d i f f • 

difusion through the transition sphere and the life time, t, , of RND. By 

life J 

adopting the renormalized vaporization laws, Eqs. (8) and (10), and a standard 
diffusion time for the transition sphere of RND, one concludes that the Q — S 
approximation is valid when the following inequality is satisfied 
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( 33 ) 


= 2C n? -2- «. (1 + B) 
t life V tS p £ n 


where B = C (T-T )/L. 

P b 


With p ~10^kg/m\ p~l kg/n\ T-T^~100°K L-10 Kcal/kg, 

•iif t diff . s 

C -0.5 Kcal/Kg-K the numerical values of — is 5 x 10 

P lif£ 

Alternatively, by using the asymptotic form of for high , deduced 
from Eq. (30), i.e. , c ~G~i = (1 G r giOT* one transforms Eq. (33) into the 

V K.N O S t S 

following form 


gpG > > G* (34) 

b s s 

where G =12 — — In (1 + B) 

S P £ 

2 

Recalling that G =4nnr„ is equal to a third of the void fraction, and 
° s y, o 

gp~0(l), one concludes that the Q-S approximation is valid in high G-sprays 

* 

when the void fraction exceeds a third of the values of G g . For typical values 

of p , p , T-T, , C and L given above, Eq. (34) reduces to 
Z b p 

nr^ > > 4x10 ^ (35) 

lo 

For example, with SMD=100y , the Q-S approximation is valid when n is greater 
than 3200 drops/cc, whereas for SMD=200y , the corresponding number density is 
400 drops/cc. 

The criterion (33) or its alternative form can be encoded in a spray 
numerical code to test the validity of Q-S assumption at each point of a spray 
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flow field. When the approximation is invalid, an alternative transient 
droplet laws should be used to determine the interfacial processes. Note that 
such droplet transient processes should be able to predict the change in the 
nature of interfacial processes. For example, a vaporizing droplet may reach 
the state of saturation when the droplet spacing falls below the critical 
value or when the thermal shielding retards the heat flux to the droplet. 
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ELEMENTAL DROPLET MODELS 
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Fig. 1 Droplet Models; Elemental Droplets 




COMPOSITE DROPLET MODELS 
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Fig. 1 Droplet Model; Composite Droplet 
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Fig. 5 


y- stratification; Vaporization 

in transition sphere for three selected test droplet 

1 _ -•! nnr 



log((T-Tb)/(Tl-Tb)) 
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Fig. 6 T- stratification = gas tenperature distribution in 
transition sphere for three selected test droplet 
bubble sizes. 
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Fig. 7 Correction factor for vaporization rate vs group com - 
bus t ion number for three selected test droplet size; 1 * 
Exact Solution. 
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